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Abstract 

The quantum Fisher information is a Riemannian metric, defined 
on the state space of a quantum system, which is symmetric and de- 
creasing under stochastic mappings. Contrary to the classical case 
such a metric is not unique. We complete the characterization, initi- 
ated by Morozova, Chentsov and Petz, of these metrics by providing 
a closed and tractable formula for the set of Morozova-Chentsov func- 
tions. In addition, we provide a continuously increasing bridge between 
the smallest and largest symmetric monotone metrics. 

1 Introduction 

In the geometric approach to classical statistics the canonical Riemannian 
metric is given by the Fisher information, and it measures the statistical dis- 
tinguishability of probability distributions. The Fisher metric is the unique 
Riemannian metric contracting under Markov morphisms 

In quantum mechanics the probability simplex is replaced by the state 
space of density matrices (positive semi-definite trace one matrices), and 
Markov morphisms are replaced by stochastic mappings. A linear map 
T: M n (C) — > M m (C) is said to be stochastic if it is completely positive and 
trace preserving. Since stochastic mappings (like their Markovian counter 
parts) represent coarse graining or randomization, one would expect statis- 
tical distinguishability of states to decrease under stochastic mappings. 

These considerations lead Chentsov and Morozova to define a monotone 
metric as a map (or rather a family of maps) p — > K p from the set Ai n 
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of positive definite n x n density matrices to sesquilinear 1 forms K p (A, B) 
defined on M n (C) satisfying: 

(i) K p (A, A) > 0, and equality holds if and only if A = 0. 

(ii) K p (A, B) = K p (B*, A*) for all p G M n and all A, Be M n (C). 

(iii) p — > K p (A, A) is continuous on M. n for every A G M n (C). 

(iv) tf T(p) (T(A),T(A)) < K„(A,A) for every p G M n , every A G M n (C) 
and every stochastic mapping T : M n (C) — > M m (C). 

It is understood that these requirements should hold for all n and m. The 
condition (ii) is sometimes omitted, but we shall only consider symmetric 
metrics. Since condition (iv) implies unitary covariance we may in all calcu- 
lations assume that p is a diagonal matrix. Chentsov and Morozova proved 
that there to each monotone metric K is a positive function c(A,p) defined 
in the first quadrant and a positive constant C such that 

n 

(1) K P (A,A) =C^Ar 1 |^r + ^|^| 2 c(A,,A,) 

for each diagonal matrix p G M. n with diagonal (Ai, . . . , A„) and every A in 
M„(C). The metric is therefore fully described by the so called Morozova- 
Chentsov function c which is symmetric in its two variables and satisfy 

c(A, A) = CA _1 and c(t\,tp) = t _1 c(A, p) 

for all t,X,p > 0. At the time it was completely unsettled which type of 
functions c would, through the formula give rise to a monotone metric. 
It was not even clear that there existed a single Morozova-Chentsov function 
(and thus a monotone metric) although Morozova and Chentsov had a few 
candidates. 

Petz connected the theory of monotone metrics with the theory of con- 
nections and means by Kubo and Ando |H] and was able to prove that the 
set of Morozova-Chentsov functions are given on the form 

(2) C(A ^ ) = _L_ A ,„ >0 , 

where / is a positive operator monotone function defined on the positive 
half-axis satisfying the functional equation 

(3) /(*) = tfir 1 ) t > o. 

1 We use the complexification proposed by Petz [12] . 
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Petz was by this result able to give some examples of Morozova-Chentsov 
functions, and they included the candidates put forward by Morozova and 
Chentsov. The existence of monotone metrics was then established. 

It is however a problem that the class of operator monotone functions 
satisfying (j3J) is largely unknown. The aim of the present paper is to provide 
a closed and tractable formula for the set of Morozova-Chentsov functions and 
in this way complete the characterization of (symmetric) monotone metrics 
given by Morozova, Chentsov and Petz. The main result is the formula given 
in Theorem 12.21 In addition, we provide a continuously increasing bridge © 
between the smallest and largest (symmetric) monotone metrics. 



2 Statement of the main results 

Let / be a positive operator monotone function defined on the positive half- 
axis. It has a canonical representation [H] of the form 

r°° f(i + c) 

(4) f(t)= / \^-d»(s) t>0, 

Jo V-r s 

where [/, is a positive (non- vanishing) finite measure on the extended half-line 
[0, oo]. The function /# defined by setting 

f#(t)=tf(t- 1 ) t G R+ 

is operator monotone. This follows easily from the above integral represen- 
tation, but may also be inferred by much simpler algebraic arguments 
2.1. Theorem (v) and 2.5 Theorem] without the use of Lowner's deep theory. 
Since = / the operation / — > /# is an involution on the set of positive 
operator monotone functions defined on the positive half-axis. The harmonic 
mean is separately (operator) increasing, hence also the function 

(5) f(t) = H(f(t),f*(t)) 2f{l)J ' U) 



f(t) + f#(t) 



is operator monotone. It is an easy calculation to show that (/)# = /. 
The formula (J^J) therefore associates a positive operator monotone function 
/ satisfying the functional equation Q to any positive operator monotone 
function /, and / = / if already /# = /. This procedure is implicitly applied 
in [121 Formula (12)] where Petz calculates a Morozova-Chentsov function 
from an operator monotone function not necessarily satisfying the functional 
equation ©. 
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There are several problems with this method, although it may be useful 
to calculate explicit examples. Firstly, the mapping / — ► f is not injective. 
There are in general infinitely many operator monotone functions which are 
mapped to the same function. Secondly, it seems difficult to specify the 
induced equivalence relation on the set of positive finite measures on the 
extended half-line [0, oo]. 

Theorem 2.1. Let /: R + — > R + be a function satisfying 

(i) f is operator monotone, 

(n) f(t) = tf\t- 1 ) for all t > 0. 
Then f admits a canonical representation 

(6) /(t)= ^ atp jf'^|. (A+ y +At) h (A)^ 

where h : [0, 1] — > [0, 1] is a measurable function and j3 is a real constant. 
Both (3 and the equivalence class containing h are uniquely determined by f. 
Any function on the given form maps the positive half-axis into itself and 
satisfy (i) and (ii). 

Proof. The result follows by applying Lemma HH] and Theorem 14.41 QED 

Note that exp/3 = f(i) exp(— z7r/4). We may adjust the constant (3 such 
that /(l) = 1. This corresponds to setting the constant C = 1 in formula 
JTJ) and gives rise to a so called Fisher adjusted metric. We are now able to 
calculate the set of Morozova-Chentsov functions. 

Theorem 2.2. A Morozova-Chentsov function c admits a canonical repre- 
sentation 

7 c(x, y = exp / — — -■- — -h(\)a\ 

K ' ^ ! x + y [ J A 2 + l (x + Xy)(Xx + y) KJ 

where h : [0, 1] — > [0, 1] is a measurable function and Co is a positive constant. 
Both Cq and the equivalence class containing h are uniquely determined by 
c. Any function c on the given form is a Morozova-Chentsov function. 

Note that c is increasing in h and that the constant Co may be adjusted 
such that c(x,x) = x' 1 . For Morozova-Chentsov functions with a fixed 
constant C and h as in formula we have c s h + (i_ s ) 9 = c^c^~ s , < s < 1. 
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Proposition 2.3. Let the exponent 7 G [0,1]. The functions 



are operator monotone, normalized in the sense that f(l) = 1 and satisfy the 
functional equation f{t) = tf(t~ r ) for t > 0. 

Proof. We first calculate the integral 

X A 2 -1 1 + t 2 „ , 2f 

CM = lOE 



/ A 2 + l (X + t)(l + Xt) to (l+t) 2 
and then by setting h(X) = 7 in © obtain the operator monotone function 



V2 Ul+() ; 

satisfying the functional equation ^j. The result now follows by setting 
= (7 -1/2) log 2. ~ QED 

Remark 2.4. Once we found the functions f 1 above, we may directly verify 
that they are operator monotone and satisfy the functional equation (JBJ) 
without making use of Theorem 12.11 Indeed, it is trivial that they satisfy 
equation To see that they are operator monotone we take a complex 
number z = rexp(i9) in the upper half plane, that is < 6 < ir. Since 1 + z 
is translated one unit to the right as compared with z it is still located in 
the upper half plane, but the angle with the real axis has decreased. It can 
therefore be written on the form 1 + z = riexp(i#x) where < Q\ < 6 and 
we notice that 0<# — 9i < 9 < n. The analytic continuation of / 7 to z may 
thus be written on the form 

/ 7 0) = 7 exp(z( 7 fl + (1 - 27)00). 

But since 76* + (1 - 27)^ = 7(6* - X ) + (1 - 7)6*1 e (0,vr), we derive that 
the imaginary part of f(z) is positive. But this proves that / 7 is operator 
monotone by Lowner's theorem [HIE]- 

A third proof is obtained by considering the set E of exponents 7 G 
[0,1] such that / 7 is operator monotone. Since /( 7 +<5)/2 = {f-yfs) 1 ^ 2 an d the 
geometric mean is operator increasing, we derive that E is mid-point convex, 
and since E is closed and contains and 1, we obtain E — [0,1]. 
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The set of operator monotone functions / defined on the positive half- 
axis such that /(l) = 1 and f(t) = for all t > has a minimal and a 
maximal element |13[ 18]. These extremal functions are given by 

f 1 ® = ^ min ^ and = (max). 

Since 4t(t + l) -2 < 1 for all t > we deduce that 

(8) < 7 < 6 < 1 / 7 (t) > f s (t) Vt > 0. 

The family (/ 7 (i)) 7 e[o,i] therefore provides a continuously decreasing bridge 
between the above extremal functions. Note also that fx/zif) — t 1 ^ 2 . The 
Morozova-Chentsov functions corresponding to the family (/ 7 (t)) 76 [ 0) i] are 
given by 

(9) c 7 (x, y) = *" V 7 (Z±A ^ 7 e [0, 1], 



and they provide a continuously increasing bridge between the smallest and 
largest symmetric monotone metrics. Finally, since c S7+ (i_ s ),5 = c 7 cj _s for 
7, 6, s G [0, 1] we realize that the mapping 7 — > c 7 is log-affine. 



3 The exponential order relation 

We introduced in an earlier paper p] the exponential ordering ^ between 
linear self-adjoint operators on a Hilbert space by setting A<B if exp A < 
exp B. It is easily verified that z< is an order relation, and since the logarithm 
is operator monotone it follows that A ^ B implies A < B. This is expressed 
by saying that the order relation -< is stronger than < . 

We also introduced and studied the set £ of real functions F : R — ► R 
which are monotone with respect to the exponential ordering. 

Proposition 3.1. The mapping $ defined by setting 
(10) $(F){t) = expF(logt) x e R 

is a bisection of £ onto the set V of positive operator monotone functions 
defined on the positive half-axis. 

Proof. Let A and B be positive invertible operators. Then 

$(F)[A) < ®(F)(B) expF(logA) < expF(log J B) 

<£> F(logA) r< ^(log-B), 

and the assertion follows since the logarithm maps the positive half-line onto 
the real line. QED 
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The next result was proved in Theorem 2.3]. 

Theorem 3.2. A non-constant function F : R — > R belongs to the class £ 
if and only if it admits an analytic continuation into the strip {z G C | < 
^sz < 7r} which leaves the strip invariant. 

Based on this result and by applying the theory of analytic functions we 
obtained Theorem 2.4] the following representation theorem. 

Theorem 3.3. A function F: R — > R is in the class £ if and only if it 
admits a canonical representation 

(11) F W=3+ £( I -i_-^) ft( A)« * £ R, 

where h : (— oo, 0] — > [0,1] «s a measurable function and (3 is a real con- 
stant. The constant j3 and the equivalence class containing h are uniquely 
determined by F. 



4 The functional equation 

The next result is the key observation in the present article. 

Lemma 4.1. Let F be a function in £. The function f = <$>(F) G V satisfies 
the functional equation f^j if and only if F(x) = x + F(—x) for every x G R. 

Proof. By taking the logarithm in the equation 

expF(logt) = $(F)(t) = /(*) = tfir 1 ) = texpFQogr 1 ), 
we realize that the functional equation Q for / is equivalent to 

F(logt) =logt + F(-logt), 

or to F(x) = x + F(—x) by setting x = logt. QED 

We need the following lemma as a preparation to the main theorem in 
this section. 

Lemma 4.2. 

2 ^ il " dx = e o<e<n. 



l A 2 -2Acos^ + l 



7 



Proof. Since the integrand can be written as 23(A — e ld ) 1 the integral is 
calculated to be 

iO iO/2 

23 [log(A - e»)] ° = 23 log — = 23 log - 

L J 1 1 + e w 2cos{9/2) 

= 23 (- log 2 cos(#/2) + log e ie l 2 ) 
= 0, 

where we used the complex logarithm. QED 

Theorem 4.3. A function F <E £ with canonical representation as given by 
l (77|) satisfies the functional equation 

F(x) = x + F(-x) \/x G R, 

if and only if h{\^ x ) = 1 — h(X) for almost all A G [—1, 0). 

Proof. Suppose that a function F G 8 satisfies the given functional equation. 
Applying analytic continuation into the strip {z G C | < 3z < n} and 
setting x = i9, we thus obtain 

(12) F(i9) = i9 + F(-i9) 0<9<n. 

Inserting the integral expression (fTT|) we then get 



° 2M h(X)dX 



CO 



A 2 -2Acosfl + l 



i9, 



or equivalently 



< 13 > L v-^ + i Wa-o o<e<«. 

We split the range of integration at the point A = — 1 and make the variable 
change A — ► A -1 in the first term and calculate 

, = C „ h(X) d\ 

1 A 2 -2Acosfl + l v ; 



CO 



-/i(A 1 ) TT dA+ / ^ — rr t — -h(X)d\. 



2A" 1 cos^ + l v A 2 y_! A 2 - 2Acos# + 1 



By applying Lemma f4. 21 we therefore obtain 
"° 2sin# 



-(h(X) + h(X- 1 )-l) 



x A 2 — 2A cos( 



which is simplified to 
1 



! A 2 - 2Acos# 



t (MA) + ma- 1 )-i) 



rfA = < < 7T. 



We introduce the change of variable u = 2A(A 2 + 1) 1 and note that u(—l) = 
— 1 and u(0) = 0. Since the derivative 

2(1 — X 2 ) 

AX) = j^j^ > - 1 < A < 0, 

we may write A = X(u) as an increasing function of u. By introducing the 
function 

0(A) = 2 a - a«) (/l(A) + h{x ~ 1] ~ 1] ~ 1 - A K °' 



we may write the above integral on the form 

1 



g(X)u'(X) dX, 



.! l-2A(A 2 + l)- 1 cos^ 

hence 

/■" I 

■^(A(u))du=0 < < 7T. 



1 

_! 1 — It COSl/ 

Setting £ = cos 9 we obtain that the function 

oo 



n=0 

hence the derivatives 



oo „ 

<^(£) =Tf / w n ^(A(«)) = 



- 1 < £ < 1, 



<p( n \0)=n\J u n g( y X(u))du = n = 0,1,2,.... 

We conclude that the function u — > g(X(u)) vanish for almost all u, and since 
A — > w(A) maps sets with positive Lebesgue measure to sets with positive 
Lebesgue measure, we derive that g(X) = for almost all A G [—1,0]. But 
this shows that /i(A _1 ) = 1 — h(X) for almost all A G [—1, 0). 

If on the other hand this relationship is assumed, we may calculate back- 
wards and obtain that F satisfies the functional equation (JT2~|) . The assertion 
then follows by applying analytic continuation and setting 6 = —ix. QED 
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Theorem 4.4. A function F : R — > R «s in the class 8 and satisfy the 
functional equation F(x) = x + F(—x) for all x G R i/ and only if it admits 
a canonical representation 

_,, . _ . l + expx Z" 1 A 2 - 1 l+exp2x , ... 

F(x) = p + log ^— + / — — - • . ./ -h{X) dX 

v ; ^2 7 A 2 + l (A + expx)(l + Aexpx) v ' 

where h : [0, 1] — > [0, 1] is a measurable function and (3 G R. TTze equivalence 
class containing h is uniquely determined by F, and [3 = 3?F(z7r/2). 



Proof Take a function F E £ with canonical representation as given by (jTljl 
and satisfying the functional equation. Then /i(A~ x ) = 1 — h(X) for almost 
all A G [—1, 0) by Theorem 14.41 By splitting the integral at the point A = — 1 
in the integral representation (|TT| and making the substitution A — ► A -1 in 
the first term, we obtain 

+ f f W-l) ft < A > dX 

J_i\X — expx X z + 1 / 

= " + / 1 (^&i-ir^) (1 -* W),tt 

+ /°(a^-a^i)^)^ 

where we used Theorem 14.31 Consequently 

+ / ( 1 \ TFTT rfA 

y_! \ 1 — A exp x X 2 + 1 / 

+ U - exp* " 1 - Aexpx " PTT + T+A* J * (A) dA ' 

and since 

A 1 I : — | dX = [ : r- — rr- d X 

_! Vl-Aexpz 1 + A 2 y J (1 + Aexpa;)(l + A 2 ) 

1 i A=1 

log(l + A exp x) — - log(l + A 2 ) 

Zi J A— 

1/1, \ X , n i l+expx 
log(l + exp ar) - - log 2 = log -j= — 
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we obtain 



. . „ 1 + expx f° 1 — A 2 1 + exp2a; , , „ 

= ^ bg + L TTa? ' (A- expI ) (1 P - A ex P ,) ft ' A > ^ 

Defining ft : [0, 1] — ► [0, 1] by setting ft, (A) = ft(— A) we may write 

n/ \ n i 1 + expx f 1 A 2 - 1 l+exp2x .... 

F(x) = ^ + log — J- + ^ ^ ■ (A + expl)(1 + Aexpl) MA) «» 

which is the desired expression. Calculating backwards we first extend h to 
the interval [—1, 0] by setting ft(— A) = ft(A), and then to the interval ] — oo, 0] 
by setting ft(A _1 ) = 1 — ft (A). We arrive in this way at the integral expression 
{TTJ), and the sufficiency thus follows by Theorem 14.31 QED 
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